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In this paper we show how the classification of topological phases in insulators and superconductors 
is changed by interactions, in the case of ID systems. We focus on the TR-invariant Majorana chain 
(BDI symmetry class). While the band classification yields an integer topological index k, it is 
known that phases characterized by values of k in the same equivalence class modulo 8 can be 
adiabatically transformed one to another by adding suitable interaction terms. Here we show that 
the eight equivalence classes are distinct and exhaustive, and provide a physical interpretation for 
the interacting invariant modulo 8. The different phases realize different Altland-Zirnbauer classes 
of the reduced density matrix for an entanglement bipartition into two half-chains. We generalize 
these results to the classification of all one dimensional gapped phases of fermionic systems with 
possible anti-unitary symmetries, utilizing the algebraic framework of central extensions. We use 
matrix product state methods to prove our results. 



INTRODUCTION 



The experimental observation of topological insulators 
in two and three dimensions has led to a renewed in- 
terest in the topological properties of insulators and su- 
perconductors. A classification encompassing all such 
band systems has been obtained by [l[ and independently 
by @- Both approaches use sophisticated mathematical 
tools, but the topological invariants they define have a 
clear underlying physical interpretation: they measure 
the twisting of the band structure over the Brillouin zone. 
The classification of and is exhaustive, but has the 
drawback that it applies only to band systems, in that 
it defines topological invariants only in a single parti- 
cle framework. In some cases, including for example the 
quantized Hall conductivity or the magneto-electric sus- 
ceptibility, these invariants can be interpreted as physical 
response functions and hence defined even in the presence 
of interactions, but in other cases there is no such inter- 
pretation. In fact, Q found a specific one dimensional 
example, the so-called Majorana chain with an unusual 
time reversal symmetry (TRS) squaring to +1, where the 
band classification is broken by interactions, in that some 
phases which are distinct in the band classification are 
actually connected in interacting Hamiltonian space. In- 
deed, while the band classification gives an integer topo- 
logical index k, it is at most only equivalence classes of 
k modulo 8 that define distinct interacting phases: Z is 
broken down to Zs. The TR-invariant Majorana chain 
example thus highlights a deficiency of the band classifi- 
cation. 

In this paper we re-examine the TR-invariant Majo- 
rana chain, and provide an interacting interpretation for 
the eight different phases. Before delving into the con- 
struction, let us first provide some intuition. The band 
Z invariant can be interpreted as the number of gapless 
Majorana modes localized at an endpoint of the chain. 
To see how it is broken, we proceed by analogy: con- 
sider an isotropic, TR-invariant spin-I chain, which has 
two phases - a trivial phase, and a topological "Haldanc" 
phase. The spin chain can also have gapless spins local- 



ized at an endpoint, but different values of this edge spin 
do not correspond to different phases: rather, it is only 
whether the edge spin is integral or half-integral that de- 
termines the phase (trivial or Haldane) of the system. 
Stated more abstractly, it is the symmetry class of the 
edge spin - "real" (T 2 = 1) or "quaternionic" (T 2 = — 1) - 
that determines the phase [|| (see also [1, @). The Majo- 
rana chain Hamiltonian respects fermionic parity (— 1) F 
in addition to T, modifying the set of possible symmetry 
classes of the 'edge spin'. In fact, there are now eight 
phases, and they are in one to one correspondence with 
k mod 8. 

In our construction we will make use of the notion of an 
entanglement spectrum: rather than looking at systems 
with physical edges, the most convenient framework for 
us will be to simulate edges with entanglement biparti- 
tions. The advantage of this approach is that it is mani- 
festly independent of the details of an edge Hamiltonian 
and does not break any extra symmetries. Indeed, the 
study of topological phases via their entanglement struc- 
ture has a rich history [H, 043 • Of particular interest to 
us is the connection between the entanglement spectrum 
and edge mode spectrum found in [j| for the fractional 
quantum Hall effect and in [if], [l2[ for band topological 
insulators and superconductors. This connection sug- 
gests that the entanglement spectrum behaves like the 
edge mode spectrum, and in particular that we should 
study the structure of representations of the generic sym- 
metries (such as T and (— 1) F ) on it. This proves to 
be a fruitful approach: already in Q it was shown that 
the trivial and Haldane phase in the spin-1 Hciscnberg 
chain are distinguished by T 2 = ±1 on the entanglement 
spectrum. Here we generalize the construction to an arbi- 
trary fermionic chain with a real time reversal symmetry, 
i.e. the TR-invariant Majorana chain. The final result is 
rather elegant: the phases of the TR-invariant Majorana 
chain are in one to one correspondence with eight of the 
ten Altland-Zirnbauer (AZ) classes fl3l - [T5f (used also in 
the scheme of [J [r| ) , and the signature of the phase one 
is in is the AZ symmetry class of the reduced density 
matrix of half of an infinite chain. 

As in the key technical tool that allows rigorous 
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arguments is that of matrix product states (MPS). The 
powerful entropy scaling bound [I?) for gapped one di- 
mensional systems allows us to approximate the ground 
state of any such infinite chain by a MPS of fixed bond 
size-, depending only on the gap and desired accuracy of 
approximation. The entanglement spectrum for a MPS 
is simply the bond Hilbert space, and known results fl8j 
classify the possible representations of the global symme- 
tries on the bond Hilbert space - they are so-called projec- 
tive representations, and for the TR-invariant Majorana 
chain they can be used to construct invariants that are 
in one to one correspondence with the Altland-Zirnbauer 
symmetry classes. Our arbitrarily accurate MPS approx- 
imations then extend this result to general gapped chains. 

The remainder of the paper is organized as follows. In 
section (JTTJ) we establish notation, define the TR-invariant 
Majorana chain, and review its band classification. In 
section (jllip we review necessary facts about matrix prod- 
uct states, and show how symmetries lead to projective 
invariants. In section (|IV|) we use these results to clas- 
sify the phases of the interacting TR-invariant Majorana 
chain. In section ([V} wc extend our scheme to the gen- 
eral classification of the phases of one dimensional gapped 
systems with both unitary and anti-unitary symmetries. 
The classification uses the algebraic notion of central ex- 
tensions, which is a precise way of defining projective in- 
variants. For completeness, in section (|VI[) we compute 
the invariants for flat-band models representative of the 
8 phases of the TR-invariant Majorana chain, showing 
that all possibilities are realized, and relate them to the 
value of the topological index k modulo 8. We conclude 
with a discussion of related matters and future directions 
in section (|VII[) . In the appendices we discuss in more 
depth the mathematical structure common to Altland- 
Zirnbauer theory and the TR-invariant Majorana chain, 
and give a short review of semisimple algebras. 

As we were writing this article, we learned of the si- 
multaneous independent work of [19( , whose results agree 
with ours. Furthermore, after the completion of this 
work, we noticed (20| , which also points to a similar clas- 
sification of gapped phases of one dimensional chains. 

II. PRELIMINARIES 

We start with a second quantized Hamiltonian H in 
the creation and annihilation operators aj , aj of spinless 
fermions, where j indexes sites of a chain. We assume H 
is gapped and includes only short range bounded strength 
interactions. We allow possible pairing terms, e.g. ajaj_. 
In this basis, the time reversal symmetry T acts as com- 
plex conjugation of the wavefunction, while the creation 
and annihilation operators are real: 

Ta 3 T- x = Oj-, 

Ta]T- x = a]. (1) 

Let us now restrict to quadratic Hamiltonians and re- 
view the band classification of the phases of the Majo- 



rana chain. It is convenient to introduce the Hcrmitian 
"Majorana" operators 

c 2 j-i = -i(aj - a]) 

C2j = ^ + a] (2) 

Any Hermitian quadratic Hamiltonian H can then be 
written in the form 

H =l ^ A; m C;C mi (3) 
l,m 

where A is a real anti-symmetric matrix. 
H : 



1 2 3 4 5 6 7 8 

Hi = 



1 234567 8 

FIG. 1: Schematic representation of the Hamiltonians Ho and 
Hi. The dots denote Majorana fermions, and the edges the 
quadratic couplings in the Hamiltonian. 

First, assume that we do not impose T: there are then 
two distinct phases of the Majorana chain Q: the trivial 
phase 

. N 

Ho = -^Y^c 2 j-ic 2j (4) 
3=1 

= E(°5°i-^) 

3 

where all sites are decoupled and unoccupied in the 
ground state, and a non-trivial phase 

. N-l 

H i = -^E c 2j c 2j+i ( 5 ) 

3=1 

= 5 E (~ a ] a J' +1 ~ a i+i ^' + a \ a ]+i + a i+ 1 a i) ■ 
3 

These are illustrated in figure [TJ The defining character- 
istic of the non-trivial phase is a two-fold ground state 
degeneracy: the edge "dangling" Majorana operators Ci 
and C2N can be paired up into a physical fermion mode 
that does not cost any energy to occupy. 

Now impose T. There are then infinitely many phases, 
indexed by an integer k 0] ■ We can interpret these phases 
in terms of dangling Majorana modes as well. We have 
infinitely many phases, instead of just two, because T 
symmetry restricts the allowed interactions, preventing 
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demonstrate via explicit calculation for Ho ((4]) and Hi 
([5]) . The Jordan- Wigner transform of Hq is 



(8) 



whose ground state of all spins pointing up is an eigen- 
state of P. The Jordan- Wigner transform of ([5]), on the 
other hand, is the Ising Hamiltonian 



(9) 



FIG. 2: Schematic representation of the Hamiltonian Hk- The 
dashed vertical line can represent a physical cut or an entan- 
glement bipartition. Note that on the left side of the cut we 
have k unpaired Majorana fermions, C2, . . . , c 2 fe. 



us from gapping out dangling Majorana modes in pairs. 
For example, consider the Hamiltonian Hk illustrated in 
figure [U and take k = 2. If we cut the chain along the 
dashed line, the unpaired Majorana modes c 2 and C4 on 
the left side cannot be gapped out, because the interac- 
tion i C2C4 is not T invariant by virtue of 2 and 4 having 
the same parity. Thus H 2 represents a new T-protected 
phase, distinct from the trivial one. One can have phases 
with any number k of such dangling Majorana modes 
whose indices have the same parity, leading to an integer 
band topological index k - the corresponding Hamiltoni- 
ans Hk are represented in diagram form in figure [2] we 
think of Hk as the Hamiltonian of k parallel Majorana 
chains, which can be viewed as a one dimensional chain 
when one chooses an appropriate ordering of the Majo- 
rana sites. It is only the equivalence class of k modulo 
8, however, that is well defined in the interacting setting, 
as we demonstrate in the rest of the paper. 

In order to define the interacting invariants, we will 
need to use matrix product states (MPS). To do this, 
we must first do a Jordan- Wigner transform to a bosonic 
spin chain. This is a general procedure that works equally 
well for interacting systems: we define 



k<j 



<7 V - = —I I 



5) Hi 1 2a W 



k<j 



(6) 



The Z2 fcrmionic parity is given by the operator 

p=n(^ 2i -ic 2J )=n^- (7) 



In general, the Jordan- Wigner transforms of the Hamilto- 
nians Hk have ground states which spontaneously break 
P if and only if k is odd. This can be explicitly seen by 
expressing the exact ground states in bosonic variables, 
resulting in so-called cluster states [2l[, but is easier to 



whose two ground states spontaneously break P. We 
stress that this symmetry breaking occurs only in the 
bosonic spin chain, and is a result of the non-local nature 
of the Jordan- Wigner transformation. P is never broken 
in the fermionic Majorana chain. In the fermionic lan- 
guage, the dimers formed by c 2 j and c 2 j-i are invariant 
under the action of P, and therefore P may be replaced 
by 



P = IC\C2N- 



(10) 



Thus, the even and the odd superpositions of the bosonic 
ground states correspond to the "empty" and the "occu- 
pied" states of the pair of dangling Majorana modes. 



III. INTERACTING INVARIANTS VIA 
MATRIX PRODUCT STATES 

Let us now review some facts about MPS which we 
will need in the following sections [3, EH, HI, ■ We as- 
sume translational invariance for convenience, though we 
believe our results to be valid even without it. An MPS 
on a one dimensional chain is a quantum state whose 
Schmidt decomposition across any cut, say, between sites 
n — 1 and n, has bounded rank a max : 



E 



n R ), a q > 0. 



(11) 



Had we performed the cut between sites n and n + 1, 
we would have obtained the same Schmidt decomposi- 
tion, but with eigenvectors | x f'" L ) and We can 
write one set of Schmidt vectors in terms of the other; 
for example 



<g> m) 



(12) 



In a more abstract language, A is a linear map of type 



U L -> Hi 



H 



spin t 



where TLl is spanned by the left 



Schmidt vectors. We refer to Hl as the bond Hilbert 
space, or entanglement Hilbert space. Thus, A m : T-Ll — > 
Hl for each m. 

To find the inherent constraints on the matrices A m , 
let 



(13) 
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and 

E* : X ^Y, A m XA l- (14) 

m 

Then the properties that the two sets of Schmidt vectors 
are orthonormal and the corresponding numbers A a are 
the same imply that 

E(l)=t, E*(A 2 )=A 2 . (15) 

In other words, 1 = (S aa ') and A 2 = (XpSppr) are a left 
and right eigenvector associated with the v = 1 eigen- 
value of the transfer matrix 

E aa ,. M ,=Y,A^{A^'y. (16) 
a/ 

If this eigenvalue is simple and there are no other eigen- 
values of magnitude 1, the MPS (or rather, the set of 
matrices defining it) is called simple, or ergodic. It is 
easy to see that all eigenvalues satisfy |^| < 1 (for exam- 
ple, because E does not increase the operator norm). 
Proceeding iteratively, we can write 

!*>=(••• A?l_A^AZl +i ■■■)]... m k . imk m k+1 ...), 

(17) 

where we implicitly sum over repeated indices. The sim- 
plicity condition mentioned above implies that the in- 
finite state |W) is pure. When we approximate ground 
states of the Jordan- Wigner transformed TR-invariant 
Majorana chain, this occurs for even topological index k. 
When k is odd, there are two degenerate ground states, 
i.e. the ground state density matrix has rank 2. It can be 
approximated with a mixed MPS whose transfer matrix 
has two eigenvalue v = 1 eigenvectors. 

Since one may feel uncomfortable with infinite states, 
let us calculate the density matrix of a finite segment. 
The equations will look more symmetric if we define the 
matrices T m by 

A m = AT m . (18) 

Then, 

pn=E^I*m)«]I. ( 19 ) 
|*pfn]>= E (r mi A-Ar m /Vi-mn). (20) 

mi ,. . . ,m n 

If the MPS data are simple, then the states j^r^i) be- 
come orthonormal in the n — > oo limit, indicating that 
the segment endpoints decouple. 

For our purposes, all nontrivial properties of MPS's 
can be derived from the following result [HI, [23| . 

Theorem. Two sets of simple MPS data, (r m , A) and 
(T' m , A') define the same physical state (i.e. P[o, n ] = p{ „i 



for all n) if and only if there is some unitary matrix U 
and number 9 such that 

= e ie U~ l T m U, A' = U^AU. (21) 

The factor e 10 is unique, and U is unique up to a phase. 

The uniqueness can be shown pretty easily. Let us sup- 
pose that eq. (|21[) has two solutions, (Ui, 0i) and (U2, 62), 
and let V = UxU^ ■ Then V commutes with A and T m 
(up to a phase), and hence with A m — AT TO . Applying 
CGI and dHJ), we find that 

E(V) = e i{ei - 02) E(t) V = e <(»i-fc)y 

so that V is an eigenvector associated with a magnitude 
1 eigenvalue of the transfer matrix. Due to the simplic- 
ity assumption, 1 is the only such eigenvalue, and the 
eigenvector is proportional to the identity matrix. For 
the proof of existence, see [lH. |23|. 

Suppose now that we have a simple MPS and a fi- 
nite global internal symmetry group G, whose elements g 
act in the physical spin space by a matrix representation 
(flWi)- The corresponding matrices T m are transformed 
into T' m = J2 n Sm-nXn- In the case of anti- unitary g the 
relevant operator is the unitary matrix (g m n) composed 
with complex conjugation, i.e. T' m = ^2 n g m -nX* n - Since 
T m and T' m represent the same physical state, the above 
theorem gives an induced action of G in the bond Hilbert 
space P.[l8|. That is, for each g there is a corresponding 
operator g, whose explicit form in some basis is 

- = f U g , if g is unitary, 

^ \ U g K, if g is anti- unitary, 

where U g is a unitary matrix and K is complex conjuga- 
tion. Such g is defined up to a phase, commutes with A, 
and satisfies 

Y,9mn{gT n g- 1 )=e i6a T m (23) 

n 

for some 9 g . Because gh and gh are associated with the 
same transformation gh of physical spins, the uniqueness 
part of the above theorem guarantees that they are equal 
up to phase, showing that the operators g form a projec- 
tive representation of G on the bond Hilbert space. 

This analysis is applicable to an arbitrary gapped spin 
chain whose Hamiltonian H is invariant under G and 
has a unique ground state. Indeed, the ground state 
can be approximated arbitrarily well with MPS [I?}, al- 
lowing us to define a projective representation of G on 
its entanglement Hilbert space as well: g : Hl — > ~Hl- 
A group element g acts on operators localized on the 
left half of the chain by AT n- cr g (X) = gXg- 1 . We 
will use this action to construct invariants which char- 
acterize the gapped phases of H. In the mathematical 
literature such invariants go by the name of extensions 
of G with U(l), and (for a finite G) are characterized 
by the group cohomology with U(l) coefficients. Infor- 
mally, they encode in an invariant way the extent to 
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which the phase ambiguities inherent in the definition of 
g are incompatible with the original group structure of 
G. But different projective representations correspond- 
ing to the same extension are compatible with each other. 
For example, the entanglement Hilbert space of a spin- 
1 chain may include sectors corresponding to different 
half-integral values of the the bond spin, e.g. 1/2 and 
3/2. The matrices T n can mix them according to the 
fusion rule [1/2] x [1] = [1/2] + [3/2]. On the other hand, 
integral and half-integral values of the bond spin cannot 
coexist without breaking the ergodicity. 

For finite G, these invariants are discrete, and using 
MPS approximations can be shown to vary continuously 
as one moves within the same gapped phase, implying 
they are in fact constant in that phase. Hence they are 
many-body invariants characterizing interacting gapped 
phases of the spin chain. We will discuss their gen- 
eral mathematical structure, for both unitary and anti- 
unitary symmetries, in section (|V| ; first, however, we 
explicitly construct them for the relatively simple case of 
the TR-invariant Majorana chain. 



IV. TR-INVARIANT MAJORANA CHAIN 

In the TR-invariant Majorana chain, the symmetry 
group G is generated by the fermionic parity P and 
time reversal T. These satisfy the commutation rela- 
tions P 2 = T 2 = 1 and PT = TP. What invariants 
can we construct for the projective action of G on the 
entanglement Hilbert space? Let us begin by consider- 
ing even topological index k, where the ground state is a 
simple MPS. We have two non-trivial operators, P and 
T, as well as their product PT. The projective form of 
the relation P 2 = 1 is P 2 = e^t for some <b; however, 
4> can be removed by a phase shift of P, so that we can 
assume P 2 = 1. T and PT also square to 1, but because 
they are anti-unitary, the corresponding projective oper- 
ators have to square to the identity up to sign: T 2 = ±1, 
(PT) 2 = ±1. These signs cannot be removed by any 
phase rotations, and together define four invariants for 
the projective representations, distinguishing four even 
k phases of the TR-invariant Majorana chain. To refor- 
mulate slightly, the phases are distinguished by whether 
T is real or quatcrnionic (column D in table |T]), and 
whether it commutes or anti-commutes with the projec- 
tive fermionic parity operator P (the invariant a = ±1). 
These invariants as functions of k will be calculated in 
section I VII Now, the reduced density matrix of the 
left half of the chain, A 2 , contains both parity-even and 
parity-odd sectors, and T commuting or anti-commuting 
with P means that T preserves or exchanges these sec- 
tors, respectively. Thus, the four different phases are 
distinguished by whether the anti-unitary symmetry T is 
real or quaternionic, and whether it is to be interpreted 
as "time reversal symmetry" T + or "particle-hole sym- 
metry" T_. This determines the orthogonal, symplectic, 
and C and D BdG Altland-Zirnbauer classes, see table |U 



(The translation to the AZ language is formalized in ap- 
pendix [A]) 

For odd topological index k one can perform a similar 
analysis, but it is complicated by the fact that there are 
now two degenerate ground states, \^>') and |^"). We 
can approximate \^>') by an MPS (r^, A'). Then because 
|*") is related to 1*'} by the Z 2 symmetry f (J7J>, a good 
MPS approximation to it is given by the MPS 



(T£,A") = (J2 P mnK, A' 
V n / 

= ((-1)X.A'), 



(24) 



where m = and m = 1 refer to the "up" and "down" 
spin states, respectively. According to this definition, 
the associated operator P' on the bond Hilbert space is 
trivial, since the relation between T' and T" is entirely 
due to the action of P on the physical spin. (Later we 
will define a nontrivial P acting in a larger space.) 

Consider now the action of time reversal symmetry. 
The states \^f') and |W) obey the clustering condition 
on correlation functions in the thermodynamic limit, and 
the action of T preserves this clustering condition. Thus 
T must either fix \*f?') and \^") individually, or exchange 
them. Let us analyze each of these cases in turn. 

If T fixes the states, then the simple MPS (r^, A') has 
a projective representation T' of T on the bond Hilbert 
space, with the two possibilities: (T') 2 = ±1. By the 
definition d24} of (r^,A") we sec that these MPS data 
are preserved if the symmetry T = PTP~ X acting on 
the state \^f") is accomplished by the operator T" = 
P'f'(P')- 1 = T". In the case that T exchanges \V) with 
I'J'"), we may apply the above argument to the symmetry 
PT, which fixed each state. 

The previous arguments arc only applicable to "unbro- 
ken" symmetries. To represent other symmetries (e.g. T 
in the second case, where it exchanges |\f' / ) and \^f") ), we 
double the bond Hilbert space and consider the nonsim- 
ple MPS that corresponds to the total (mixed) ground 
state. It is the direct sum, 



r, 



r 

rn 

o 



o 

r" 



A 



A' 
A" 



(25) 



with the symmetry P acting as 



1 



P = a* ® 1 = ^ q J ■ (26) 

Furthermore, we can define an additional operator 

Z = a* ® 1 = (I , (27) 

which distinguishes the sectors of the Hilbert space cor- 
responding to |*'), |*"). We have ZP = -PZ. In this 
representation, the symmetry T acts as 



T = 







rp ( 

or T = ~, 


f'\ 


Vo 


V 


\T 


o 



(28) 
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if f fixes or exchanges the two states, respectively. We 
see that PZ = aZP, where a is equal to 1 in the first 
case and to — f in the second case. 

To translate to the AZ language, we define the "time- 
reversal symmetry" by T + = T (because T always com- 
mutes with P). The "sublattice symmetry" Z should 
be equal to T+T_, therefore T_ = f^Z. Thus we have 
obtained the chiral orthogonal and chiral symplcctic sym- 
metry classes for a = 1 , and the CI and Dili BdG classes 
for a = — 1 (see table UJ. 

V. GENERAL CLASSIFICATION 

In this section we consider an arbitrary gapped one 
dimensional system that consists of spins or itinerant 
fcrmions and has an unbroken symmetry group involving 
both unitary and anti-unitary elements. (Let us reiterate 
that the fermionic parity P is also always unbroken - the 
Z2 breaking discussed in the previous sections was in the 
Jordan- Wigncr transformed bosonic spin chain, which is 
related in a non-local way to the fermionic chain.) The 
TR-invariant Majorana chain is in some sense the sim- 
plest non-trivial example, in that the Hamiltonian is in- 
variant under both the unitary fermionic parity P and 
anti-unitary time reversal T. Having shown how to rig- 
orously define the invariants for the TR-invariant Majo- 
rana chain with entanglement bipartitions, we shift our 
point of view slightly in this section, and discuss systems 
with physical edges - that is, finite segments. We will be 
concerned with low energy degrees of freedom localized at 
the endpoints of the segments. Although all the proper- 
ties we need can be proven by translating the statements 
to the language of entanglement bipartitions and using 
MPS as in the previous sections, here we will be con- 
tent with simply clearly stating the necessary properties. 
These will have to do with the form of the low energy 
sub-spaces, and the space of operators acting on them. 

We begin by considering a bosonic spin system defined 
on a one dimensional lattice, with a Hamiltonian H that 
is gapped in the bulk and is invariant under the simul- 
taneous action of a symmetry group G on all sites. This 
action may be projective, which allows us, for example, 
to treat a (dimerized) spin-1/2 Heisenbcrg chain as a sys- 
tem with SO(3) symmetry In any case, there is a genuine 
action of G on physical obscrvables, and that is sufficient 
for our purposes. Now we imagine a finite chain, with left 
and right endpoints. The Hilbert space £bound of low en- 
ergy boundary states is then assumed to decompose as 

Abound = Ci® C r . (29) 

We find it more useful to consider the algebras of linear 
operators defined on these Hilbert spaces: 

C(£ bound ) =C(Ci) ®C(Cr). (30) 

The group G has a well defined action by automorphisms 
on such operators. 



Let us specialize to the right endpoint, and denote the 
action of g 6 G on A r = C(C r ) by X >-> <j g {X). We first 
consider only unitary symmetries. The map a g being an 
automorphism, preserves scalars and satisfies 

<J g (X + Y) = a g (X)+a g (Y), 
a g (XY) = a g (X)a g (Y), 
a g (Xi) = <j 9 (X)K (31) 

Because A r is a matrix algebra, and hence simple, every 
automorphism must be of the form cr g {X) = U g X U g x 
for some unitary U g , which is well defined up to a phase. 
To accommodate future generalizations, let us use an al- 
ternative notation: g — U g . We fix the phase for each g 
in an arbitrary way, but also consider an invariant object, 
the group G that consists of operators of the form e 1 ^ g. 
Thus we have an exact sequence 

1 ->• U(l) ->• G ->• G ->• 1. (32) 

The isomorphism classes of such exact sequences are 
known as extensions of G with U(l). For example, the 
group G = SO(3) has two extensions: G = SO(3) x U(l) 
or G = U(2). Representations of the two variants of 
G correspond to integral and half-integral spin, respec- 
tively, where we require that the U(l) subgroup act by 
the multiplication by scalars. The first case is realized 
when the spin-1/2 chain is cut between dimers, whereas 
the second case describes the chain cut across a dimcr. 
Another interpretation is the trivial and the Haldane 
phase in a spin-1 chain. As shown in Refs. [3, [H|, the 
two phases remain distinct when the symmetry is re- 
duced to the dihedral subgroup D 2 C SO (3). Indeed, 
£>2 = Z2 x ^2 = {/, X, y, XY} also has two extensions. 
The nontrivial one is given by, e.g., Ui = 1, Ux = <? x , 
U Y = cry, Uxy = cr z . 

In the following discussion we take G to be finite. Ex- 
tensions of G with U(l) can be described explicitly using 
2-cochains, i.e., U(l)-valued functions C(g, h). These are 
defined by 

gh = C{g,h)g~h. (33) 

Associativity gives the cocycle constraint, 

C(g, h) C(fg, hy 1 C(f, gh) C(f, g)- 1 = 1, (34) 

while redefinition of the U g by a U(l) phase leads to the 
gauge symmetry C ~ A • C, where the right-hand side is 
the product of C(g, h) and A(g, h) = B(h) B(gh)" 1 B{g) 
for some function (1-cochain) B. Gauge equivalence 
classes of functions C(g, h) satisfying the associativity 
constraint constitute the cohomology group H 2 (G, U(l)). 
As we have seen, they are in one to one correspondence 
with U(l) extensions of G. 

With a slight modification to this framework we can 
incorporate anti-unitary symmetries. The setup now is 
that G is equipped with a homomorphism a : G — > Z2, 
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indicating whether a given element g is unitary or anti- 
unitary. The action of G on scalars is generally nontriv- 
ial, 



09(c) 



if a(g) 
if a(g) 



-1, 
-1, 



(35) 



while equations (J3Tj) still hold. Each automorphism a g is 
still represented by an operator g acting in the space of 
boundary states C r . If we fix a basis, we can write g = U g 
or g = U g K , where U g is a unitary matrix and K is com- 
plex conjugation. The gapped phases are still classified 
by U(l) extensions, or elements of H 2 (G, U(l)), where 
the action (f35|) of G on U(l) is assumed. This amounts 
to a simple modification of the cocycle condition, 



a f (C(g, h)) C(fg, ft)" 1 C(f, gh) C(f, <?)" 



(36) 



and the equivalence relation: C ~ A ■ C , where 

A(g,h) = a g (B(h))B(gh)- 1 B(g). (37) 

The last equation corresponds to the gauge transforma- 
tion g 1— > B(g) g (the multiplication order matters). One 
example, G = {I, T} with T anti-unitary, was mentioned 
in the introduction. 

The incorporation of fermions is slightly more involved. 
We now have a special unitary symmetry, the fermionic 
parity P £ G which is involutory (P 2 ~ I) and central 
in G. Operator algebras are Z 2 -graded, ^ = ^(°)©^( 1 ). 
Namely, an operator X is even (odd) if it preserves (resp., 
changes) the number of fermions modulo 2, i.e., 

XeA (x) iff <Jp{X) = {-l) x X (a; = 0,1). (38) 

At the elementary level, such operators are sums of even 
(odd) products of aj , . Note that all automorphisms 
a g respect the grading. Indeed, if X £ A^ x \ then 

a P (a g (X)) = a g (a P (X)) = (-l) x a g (X), 

implying that a g (X) £ A (x \ 

The algebras Ai , A r of operators acting at the ends of 
the fermionic chain are simple in the graded sense. In- 
stead of defining this notion abstractly, we will describe 
the structure of such algebras and give some physical in- 
terpretation (for the general theory of simple Z2-graded 
algebras over an arbitrary field, see Ref. [24|). Note that 
A r consists of operators that act near the right endpoint 
in the actual system rather then its Jordan- Wigner trans- 
formed version. The MPS interpretation is slightly differ- 
ent; let us only mention that the operator Z (see eq. (j2"T|) ) 
belongs to A r , but P artificially extends that algebra. 

There are two alternatives. An even simple algebra 
(over the field of complex numbers) has the form A = 
C(£W®£W), where the even and odd parts are defined 
as follows: 



Aoo 
. 





An 



A®, 




A10 



A01 




A {1) . (39) 



'2m+1 1 



4m+3 4m+4 



4m+5 4m+6 



4m+1 4m+2 



FIG. 3: Schematic representation of the Hamiltonian H2 m +i- 



Thus, the grading automorphism up is represented by an 
involutory operator P: 



a P {X) = PXP, 



P 









-1 C (1) 



(40) 



The spaces C^°\ consist of states with different lo- 
cal fermionic parity, though we cannot tell which is which 
since P is defined up to a sign. An odd simple algebra 
does not have an internal grading operator, but rather an 
odd involutory central element Z. For example, consider 
the model 7?2m+i m figure [3] and think of the algebra 
Ai of operators that act on the degenerate ground state 
and can be localized near the left end of the chain. This 
algebra is generated by k = 2m + 1 unpaired Majorana 
modes, ci, C3, . . . , C4 m +i and has an odd involutory cen- 
tral element Zi = i m c\c^ ■ ■ -Gim+i. In general, an odd 
simple algebra has the form A = C(£) © <C(C)Z. Such 
algebras have a convenient representation on the Hilbert 
space C 2 ©£: one sets Z to M%\c, where M is a suitable 
2x2 matrix. In the previous section we had M = cr 2 , 
but it may be more natural to use a v and a x at the left 
and the right end of the chain, respectively. Thus, there 
are two standard forms of general algebra elements: 



A',, 






x 





iXi 



-iX 1 




^0 
Y 



£A(°\ 







£A^ 



A^ 1 



(41) 



(42) 



For a system composed of pairs of Majorana modes, 
the left and right algebras are either both even or both 
odd. The full algebra Abound is the graded tensor product 
Ai ®gr*4r, which implies a special multiplication rule. It 
is sufficient to define the multiplication for elements of 
the form X ® gr 1 and 1 (x) gr Y, which we identify with X 



and Y, respectively. 
0, 1), then 



EX £ A 



(x) 



and Y £ A r v) 



(x,y = 



XY = X 



Y. 



YX = {-l) xv X® gr Y. (43) 
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In other words, operators at the two ends commute or 
anticommute according to their parity. This rule is natu- 
ral if we think about physical fcrmions, but one can also 
construct such a product by modifying the usual tensor 
product, a trick appropriate in the Jordan- Wigner pic- 
ture. In the even case, all operators act in the Hilbert 
space (C\ 0) ©£, (1) ) © (£r 0) ®CP) with the obvious grad- 
ing, and the (anti)commutation is achieved by defining 
1 © gr Y as [Pi) V <8> Y. In the odd case, the Hilbert space 

is Abound = C 2 ® Ci © L r with P = a z © 1 © 1, and 
the correct commutation relations are due to the use of 
representations ([31]), (|32]l. 

Finally, we consider the action of symmetries on the 
simple Z 2 -graded algebra A = A r . In the even case, 
A = C(jC^ © C^) remains simple even if we ignore the 
grading. Therefore we can represent every automorphism 
a g by a unitary or anti-unitary operator g (in particular 
up is represented by P). This structure is just an ex- 
tension of G with U(l) in the usual sense, but in the 
fermionic context we call it an even extension. It auto- 
matically defines a homomorphism /i : G — > Z2, which 
tells whether a given symmetry g preserves or changes 
the local fermionic parity: 

geA^\ M (.g)e{0,l}. (44) 

To see that g is indeed homogeneous, we write this con- 
dition as follows: PgP = (—l)^ 9 '^, or 

a g (P) = (-ir^P. (45) 

The algebra elements P and <J g (P) are related by a sign 
because they are both involutions and represent the same 
automorphism, as shown by this calculation: 

a g (P)Xa g (P) = a g (Pa g - l (X)P) 

= a g (a P (a g -i{X))) = a P (X). 

It is clear that fi(P) = 0. 

In the odd case, the algebra A = C(£) ©C(£)P is not 
simple without grading, but the subalgcbra A^ = C(£) 
is. Therefore the restriction of the automorphism a g to 
can be represented by an operator g' G A^ . Since 
up acts on A^ trivially, we obtain an extension of the 
quotient group G' = G/{I,P} with U(l). To find the 
action of a g on the whole algebra, we also need o~ g (Z). 
It is defined by analogy with equation (|45[) , namely 

a g (Z) = (-ir^Z. (46) 

(This definition is correct because o~ g preserves the prop- 
erties of Z to be involutory, odd, and central, and such 
an operator is unique up to sign.) In the present case, 
the homomorphism (i is an independent piece of data sat- 
isfying the only constraint /Lt(P) = 1. For an arbitrary 
X e A x we have 



where g' is constant on each coset {g,gP}. 

One may wonder why the structure of even and odd 
extensions is so different. We cannot give an exact an- 
swer but notice that odd extensions do not exist for all 
groups. Indeed, a homomorphism fi : G — > Z2 with the 
property n(P) = 1 defines an embedding r : G' —> G. 
(In each coset {g, gP} we pick the element r(g) = h such 
that n(h) = 0.) Thus, G = G" x {I,P}. In particu- 
lar, P cannot have a square root in G. This precludes 
odd extensions for the groups G = 1 4m and their limiting 
case U(l), which is consistent with the fact that unpaired 
Majorana modes only occur in superconducting systems. 

On the other hand, if G = G' x {/, P}, then even and 
odd extensions can be unified. They are both defined by 
an (ordinary) extension of G' with U(l) and a homomor- 
phism /i : G — > 1*2- The latter can be specified by the 
value of n{P) and a homomorphism fi' : G' — > Z2. For 
the TR-invariant Majorana chain, we have G' = {I,T} 
where T is anti-unitary. Thus, a general extension is 
given by the invariants e = (—l)^ p \ a = (— I) M ' T ', and 
T 2 = ±1. The latter may be interpreted as a choice be- 
tween real numbers (D = R) and quaternions (D = H). 
These are the Wall invariants in table |U 



VI. CALCULATION OF THE INVARIANTS 

We have defined interacting invariants which distin- 
guish eight possible phases of the Majorana chain. We 
now show, for completeness, that all eight phases are 
realized in the band TR-invariant Majorana chain, and 
relate the band topological index k to the values of the 
interacting invariants from the previous section. Here we 
essentially re-derive the mathematical structure of the 
previous section, for the specific case of the TR-invariant 
Majorana chain, in a much more pedestrian fashion. 

We consider the flat band Majorana chain Hamilto- 
nians Hk defined in figure [5] We need to compute P 
and T on Hl , the Hilbert space spanned by the Schmidt 
cigenstates of the left half of the chain. (It is similar to 
the space of low energy states C r that was used in the 
previous section.) We begin by taking k even. 

A. Even k 

Let k = 2m. The picture is as in figure [2] where we 
take the dashed vertical line to represent the right end- 
point of the system. Then the space of low energy opera- 
tors <C(T-Ll) is spanned by the unpaired Majorana modes 
C2, C4, . . . , C2fe. First, consider the action of the parity op- 
erator P. defined as the product of physical fcrmion occu- 
pation numbers over all sites. Because the bulk Majorana 
modes are paired into unoccupied physical fermions, P 
reduces at low energy to the parity operator for the 2m 
Majoranas C2, . . . , C2k- 



<T g {X) = {-iy^*g'X{g')-\ (47) 



P = (ic 2 c 4 ) ■ ■ ■ (ic 4m _2C 4m ). 



(48) 
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Second, consider the action of T. By definition, 
Tc2iT~ 1 = c 2 i for I G l,...,2m, and TiT^ 1 = —i. 
Thus from (|48|) we see that 

fPf - 1 = aP, a = (-l) m . (49) 

Now we need to determine the sign of T 2 . In general, 
for a 2 J dimensional Hilbcrt space H' of the Majoranas 
Ci , . . . , c' 2 j , the action of an anti-unitary T 1 is uniquely 
determined up to phase by its adjoint action on the c[. 
Indeed, given two anti-unitaries T' and T" with the same 
adjoint action, the unitary (T')~ 1 T" commutes with the 
entire Clifford algebra, so that T" is equal to T" up to a 
phase. In particular, because T is anti-unitary, the sign 
of (T') 2 is uniquely determined. Now, we know that for 
even indices 21, 

TcuT- 1 = c 2l (50) 

for I G 1, . . . , 2m. Let us find an explicit expression for 
T in the occupation number basis, assuming that the 
annihilation and creation operators b ri b\ (r = 1, . . . , m) 
arc defined as follows: 

C4r-2 = -i(b r - b\) 

Cir = b r + b\.. (51) 

Let |0) be the state annihilated by all the b r . Any \tp) G 
Hl can be written as 

W = £ C ai ... am b\ ai ---bl am \0)- (52) 

a,e{0,l} 

Let K be the complex conjugation operator on this space 
of wavefunctions: that is, K acts by complex conjugating 
the coefficients C Ql ... Qm . Now define 

f _ f c 2 c 6 • • ■ c Am _ 2 K for m even, , . 

1 C4C8 • • ■ C4 m K for m odd. ^ ' 

We check that T, so defined, satisfies (|50[) . and hence is 
the unique choice, up to sign, of the projective action 
of T on "Hl. We now find by direct computation that 
T 2 = 1 for m = 0, 1 modulo 4, and f 2 = -1 for m = 2, 3 
modulo 4, so that f 2 = 1 for k = 0, 2, and T 2 = -1 for 
fc = 4,6. 

Thus we see that the sign of the commutation relation 
between P and T, and the sign of T 2 , together uniquely 
determine even k modulo 8. We gather these results in 
table H 

B. Odd k 

Now let k = 2m + 1. We proceed as in the previous 
subsection, but the situation is now complicated by the 
two degenerate ground states. We need to construct the 
Hilbcrt space Hl spanned by the Schmidt eigenstates of 
both of these ground states. The setup again is as in 



figure [2j this time with an odd number of Majoranas, 
C2, • ■ • , C4 m +2, lacking partners. This situation may seem 
problematic, as we need an even number of Majoranas 
to define a physical Hilbert space. It is remedied by the 
the existence of an extra Majorana mode Coo at infinity, 
which appears whenever we try to cut the system off at 
a finite size. Thus is the 2 m+1 dimensional Hilbert 
space of the Majoranas C2, . . . , C4 m +2, Coo- 

As before, the parity operator is just the product 

P = (ic 2 c 4 ) ■ ■ ■ (ia ){ic im +2C 00 ). (54) 

It is also useful to introduce an operator Z that is similar 
to P but does not involve c^: 

Z = i m c 2 ■ ■ •c 4m + 2 - (55) 

It squares to 1, commutes with all operators acting lo- 
cally (i.e. with combinations of C2, . . . , C4, n +2), and anti- 
commutcs with P. Under the time reversal symmetry, Z 
transforms like P in the even case: 

TZf- 1 =aZ, a=(-l) m . (56) 

However, we now have an ambiguity, because we have 
not specified the transformation of Coo under T. In fact, 
there are two possible choices for a consistent anti-unitary 
symmetry, TcqoT" 1 = ±Coo. By analogy with the MPS 
construction in section IIV1 we choose the sign so that T 
and P commute, namely 

fcoof- 1 = (-ir +1 Coo . (57) 

Similarly to (|53[) , we now pair up the Cj into annihilation 
and creation operators as indicated in the expression (|54[) 
for P, and in this basis define 

j> _ / C 2 C 6 ' • • c 4 m+2-^ for m odd , 58 ^ 
\ c 2 c 6 • • • c im+ 2CocK for m even. 

(In the second case, one can represent T as T'P, where 
T 1 does not involve Cqo . Although T" does not implement 
the time reversal symmetry on the Majorana operators, it 
acts correctly on physical observables, i.e. even products 
of c 2 , . . . ,C 4m+ 2)- 

A simple calculation shows that T 2 = (T') 2 = —1 for 
m = 1,2 modulo 4, and f 2 = (T 1 ) 2 = 1 for m = 0, 3 
modulo 4. Thus, f 2 = 1 for k = 1, 7 and f 2 = -1 for 
k = 3, 5. The results are gathered in table fl] 

VII. DISCUSSION 

We have constructed interacting physical invariants 
which distinguish among the eight T protected phases 
of the TR-invariant Majorana chain, and related them 
to Altland-Zirnbaucr symmetry classes of matrices. This 
TR-invariant Majorana chain is in some sense the 
simplest one dimensional system incorporating both 
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topological index 
k mod 8 


Wall invariants 


f2 


Altland-Zirnbauer classification 


f 




a 


D 




T- 


Z 2 


Cartan label 





R 


+ 


+ 1 


R 


+ 1 


+ 1 






AI (orthogonal) 


1 


R 


— 


+ 1 


R 


+ 1 


+ 1 


+ 1 


1 


BDI (chiral orthogonal) 


2 


R 


+ 


-1 


R 


+ 1 




+ 1 




D (BdG) 


3 


E 


- 


-1 


H 


-1 


-1 


+ 1 


1 


Dili (BdG) 


4 


E 


+ 


+ 1 


H 


-1 


-1 






All (symplectic) 


5 


E 




+ 1 


H 


-1 


-1 


-1 


1 


CII (chiral symplectic) 


6 


E 


+ 


-1 


H 


-1 




-1 




C (BdG) 


7 


E 




-1 


R 


+ 1 


+ 1 


-1 


1 


CI (BdG) 




C 


+ 


+ 1 


C 










A (unitary) 




C 




+ 1 


C 








1 


AIII (chiral unitary) 



TABLE I: The correspondence between the topological index k of a TR-invariant Majorana chain, the Wall invariants of simple 
Z2-graded algebras [2J|, and the Altland-Zirnbauer classes of free-fermion Hamiltonians [L§-[l5| (using the scheme P,[l^). Both 
classifications are explained in appendix [S] and the invariants a and D for the Majorana chain are calculated in section fVTl 



fermions and anti-unitary symmetries. Its solution led 
us to a general framework for classifying gapped phases 
of such chains with an arbitrary symmetry group G, that 
of extensions of G with U(l). 

For the case of the TR-invariant Majorana chain, one 
can tell which phase one is in by looking at the ground 
state reduced density density matrix p of a semi-infinite 
chain. The key is to study the symmetries of p considered 
as an operator on the entanglement Hilbert space, which 
is defined as the span of the, say, left Schmidt eigenstates: 
p : Hl — > T~Ll- T~L carries additional structure coming 
from a projective representation P, T of the underlying 
symmetries. In addition, when one has unpaired Majo- 
rana modes there is a symmetry Z which distinguishes 
the ground states. This structure amounts to a speci- 
fication of Altland-Zirnbauer symmetry class for p; one 
aspect of this interpretation is that the anti-unitary sym- 
metry exchanging different fermionic parity (P) sectors 
is thought of as a particle-hole symmetry (T_). 

The parity of k has a familiar physical interpretation: 
it describes the presence or absence of unpaired Majo- 
rana modes at the endpoints of a sample. What about 
the invariants which determine k modulo 4 and 8? In 
certain limits, these invariants reduce to other, more well 
known ones. For example, consider restricting to insula- 
tors, by disallowing pairing terms in the Hamiltonian. 
This eliminates the four superconducting phases with 
odd k, but still leaves four distinct phases correspond- 
ing to even k = 2m. Then even and odd m correspond 
respectively to the two phases of polyacetelene. Indeed, 
it is well known that a domain wall between these two 
phases supports a gapless mode, and carries fractional 
charge ±1/2 depending on the occupation of this mode 
[25) . Now, we can view P and T as effectively acting 
on this gapless mode, and when they anti-commute, i.e. 
when m is odd, there is a Kramers pair of edge states 
with opposite fermionic parity. In the simplest case, by 
symmetry, they must have charge ±1/2. Hence the in- 
variant that determines even k modulo 4 recovers the 



well-known physics of the Su-Schrieffer model. This in- 
variant also has a field theory interpretation, as the pres- 
ence or absence of a non-trivial 9 = tt theta term in the 
effective U(l) gauge action for the theory. 

What about k modulo 8? In trying to distinguish be- 
tween the phase of Hamiltonian H4 and the trivial phase, 
we can impose the £0(4) symmetry, which rotates the 
four chains into one another. The non-trivial phase of 
H4 has edge modes that transform as (1/2, 0) © (0, 1/2) 
under so(4) = su(2) © su(2). The generic Hamiltonian 
of this phase will split the degeneracy between these spin 
1/2's, but will retain at least one gapless spin 1/2. Hence 
the non-trivial and trivial phases are distinguished by the 
presence or absence of a spin 1 /2 edge mode. This is rem- 
iniscent of the AKLT state, and indeed, in the space of 
su(2) C so(4) symmetric Hamiltonians the non-trivial 
phase represented by H4 is precisely the Haldane phase. 

Thus the invariant k modulo 8 ties together and gener- 
alizes several well known one dimensional physical invari- 
ants. It would be interesting to generalize these methods 
to higher dimensions. 
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Appendix A: Altland-Zirnbauer and Wall classes 

In this appendix, we reveal the exact mathematical 
structure that is common to the TR-invariant Majorana 
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chain and the Altland-Zirnbauer theory. The AZ classi- 
fication is concerned with a different set of physical sys- 
tems, free-fcrmion Hamiltonians (in dimension 0). In this 
setting, the symmetry group G acts linearly on a,j, at, or 
cquivalently, on Cj. As is usual, the group elements are 
marked as "unitary" or "anti-unitary" (using a homo- 
morphism a : G — > Z2). Because both types of sym- 
metries preserve Hermicity, we have a linear action of 
G in the mode space A4 = M. 2N that consists of opera- 
tors of the form X = J2i x i c i i x l ^ The Majorana 
modes q form an orthonormal basis of M. , and the inner 
product between two arbitrary elements is given by the 
anticommutator (more exactly, \{X, Y}). Thus, each 
symmetry g £ G is represented by a real orthogonal ma- 
trix S g £ 0(2N). We are interested in characterizing 
all quadratic Hamiltonians (|3|) that are invariant under 
this action. This invariance translates to the following 
condition on the matrix A in ([3]): 



A, if g is unitary, 
~A 1 if g is anti-unitary. 



(Al) 



This problem was formulated and solved in |15| using 
the Nambu space language. (A transparent, albeit brief 
exposition can be found in [la].) The Nambu space con- 
sists of complex linear combinations of aj , a ' , or equiv- 
alently, of q. Therefore the Nambu space Af = C 2N is 
just the complexified mode space. Conversely, the mode 
space M C M is the real subspace fixed by the conju- 
gation X ^ X\ which is an anti-unitary transforma- 
tion of the Nambu space. In the conventional basis of 
aj,aj, it is written as C — UcK. To extend the invari- 
ance condition to the Nambu space one could simply use 
the same equation (|A1[) . but tradition requires that the 
equations look similar to those in many-body quantum 
mechanics. So let us introduce a "single-particle Hamil- 
tonian" iJ s .p. = jA (not bothering with basis changes) 
and redefine the symmetries, so that for anti-unitary g, 
the corresponding operator on the Nambu space is like- 
wise anti-unitary: 
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S g , if g is unitary, 
S g C ', if g is anti-unitary. 



(A2) 



Note that the g commute with C . Now, the Hamiltonian 
satisfies these conditions: 

gH s . p , g^ 1 = H s , p . for all g £ G, 



C Hs.p. G 1 



-H s . p .. (A3) 



The solution turns out to be rather insensitive to the 
symmetry group G. Indeed, suppose first that G consists 
of only unitary symmetries. The action of G then defines 
a block decomposition, 



(A4) 



where A indexes the irreps C\, and the Hamiltonian acts 
independently in each T-L\ (so that dimH> is the mul- 
tiplicity of C\). The conjugation discriminates between 
real irreps, quaternionic irreps, and conjugate pairs of 
complex irreps; each choice corresponds to a certain form 
of possible Hamiltonians for the given block. For exam- 
ple, suppose that there is only one block C ® % of quater- 
nionic type, meaning that C 2 C = — 1. But the conjuga- 
tion on the total space satisfies the condition C 2 = 1 and 
factors as C = Cc <8> C-u , hence = — 1 . The anticom- 
mutation between C-u and the Hamiltonian implies that 
H s . p . £ sp(n) (the Cartan class C). Similarly, real blocks 
and conjugate pairs of complex blocks yield Hamiltonians 
in the D and A Cartan classes, respectively. 

When physical anti-unitary symmetries are allowed, 
the number of cases increases to 10. The set of possi- 
ble Hamiltonians (within a single block) is defined by 
commutation relations with a few operators that remain 
after factoring out unitary symmetries. Such operators 
may be even or odd, depending on whether they com- 
mute or anticommute with the Hamiltonian. Up to three 
special operators may be present, which are labeled as 
follows: 

• f+: even anti-unitary, f 2 = ±1 (TRS); 

• f_: odd anti-unitary, T 2 = ±1 (PHS); 

• Z: odd unitary, Z 2 = 1 (SLS). 

(The abbreviations stand for "time-reversal symmetry" , 
"particle-hole symmetry", and "sublattice symmetry", 
though they are not concrete elements of the symmetry 
group.) Furthermore, if two of these operators exist, the 
third is defined by the identity 



Z = T+T- 



(A5) 



All 10 cases are listed in tablelH Note that in this scheme, 
the special operator C mixes with physical symmetries 
and may disappear in the unitary elimination process. 

A mathematically more appealing way to arrive at the 
same result is to work with real coefficients 1 and use 
equation (|A1|) directly. The technical trick, 

real space f-> complex space with a conjugation 

may be used at some point, but the conjugation will al- 
ways be identified as such. The parity of operators is 
defined according to their commutation with A; in par- 
ticular, the S g are even for unitary g and odd for anti- 
unitary g. We now replace the symmetry group by the 



1 Here we follow Dyson's philosophy [26l |. In fact, he also found ten 
symmetry classes, though a coarser division into only 3 classes 
appears in the final answer. This is because Dyson considers 
the action of symmetries on the Hilbert space of a many-body 
system. Unlike the Nambu space, the Hilbert space of quantum 
states has no canonical conjugation, since the multiplication by 
a phase is an intrinsic symmetry. 
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algebra A = A { ^ © A^\ where _4. (0) and A^ consist 
of formal linear combinations of unitary and anti-unitary 
elements g £ G, respectively. We call A^ ' and A^ the 
even and odd parts of A. Thus, A is a semisimple as- 
sociative Z2-graded algebra over the field of reals, and 
its action on M. encodes all relevant information about 
the symmetries. (The simplicity and semisimplicity of 
graded algebras are defined in terms of representations 
as in the ungraded case, see appendix [B] However, inter- 
twiners can be even or odd] the odd ones commute with 
_4 (0) but anticommute with _4 (1) .) 

The solution to this more abstract problem also uses 
a block decomposition. Although the blocks of the mode 
space do not factor as C\ ® T-L\, they carry two sim- 
ple graded algebras: the symmetry algebra ,4a (a simple 
component of A) and the algebra B\ of operators that 
(anti)commute with these symmetries. The "Hamilto- 
nian" for each block belongs to the odd part of B\. 
This method will be presented elsewhere. It is relatively 
straightforward, thanks to Wall's complete theory of sim- 
ple ^-graded algebras (over an arbitrary field) [241 ] . 

We now give a short summary of Wall's invariants in 
the real case. We first reproduce the original definition 
(mainly, for completeness). Then we describe it in terms 
of complex algebras with conjugation, the language most 
natural in the Majorana chain context. (See section [V] 
for the classification of complex Z2-graded algebras with 
general symmetries.) Finally, we formulate a variant that 
is most closely related to the AZ scheme, see equations 

Let A = A^ ffi A^ 1 ' be a simple Z 2 -graded algebra 
over the reals. Four invariants are defined successively. 

1. Graded center 6 = Z(A) (1 A^ , where Z(A) is the 
center of A in the usual sense. Ms a field, i.e. K 
or C. It is convenient to consider A as a graded 
algebra over {. The invariants a and D (see below) 
are trivial for complex algebras. 

2. A label e, cither + or — . In both cases, there is 
an ungraded simple algebra V with center t and an 
invertible element u G A, though they have some- 
what different properties. 

(+) ( even case): T> = A (as an ungraded algebra); 
u belongs to A^ and defines the grading, i.e. 
uXu- 1 = (-If A for any X € A {x) (x = 
0,1). 

(— ) (odd case): T> ~ A^°\ whereas u belongs to 
A^ and commutes with all X e A. It follows 
that A (1) = uV. 

3. A nonzero number a = v? S t. (It is clear that 
u 2 belongs to A^ and commutes with everything, 
hence it is in the graded center.) Because a is de- 
fined up to multiplication by a square, we may as- 
sume that a = ±1 if £ = ffi and a = 1 if t = C. 

4. The type D of the ungraded simple algebra V: real 
or quaternionic if t = R, complex if { = C. 



We can readily see that there arc 10 possibilities: 8 with 
8 = K and 2 with £ = C (see table [J). 

Let us now specialize to % = M. (the case relevant to the 
Majorana chain) and represent the above structure by 
commutation relations between some operators. To this 
end, we replace A with its complexified version Ac (by 
adjoining the imaginary unit). One can prove that Ac is 
simple as a Z2-graded algebra. Multiplying the special 
clement u by a suitable complex number, we obtain W £ 
Ac such that W 2 = 1. In the main part of the paper, we 
denoted W by either P or Z: 



W 



P in the even case, 
Z in the odd case. 



(A6) 



The invariant a is defined by these equations: 

r (W) = aW, (A7) 



W 2 = 1, 



where o~t is the complex conjugation on Ac- 

In the even case, Ac is isomorphic to C(£) as an un- 
graded algebra, and o~t can be represented by an anti- 
unitary operator T acting in C. In the odd case, this docs 
not work directly because Ac — C(£)©C(£) Z is not sim- 
ple as an ungraded algebra. To achieve the desired repre- 
sentation, we double the space £, set W = Z = a z $5> l,c, 
and introduce an auxiliary operator P = a x <E> tc- (This 
procedure is described more concretely in section II VI ) 
Now we can represent or by T = 1 <g> T' if a = 1, or 



by T 



V if 



T. The whole construction is 



characterized by the following equations. 



Even case: P 



2 



1. TP = aPT: 



Odd case: 



Z l = 1, 



TZ = aZT, 



1. 



PT = TP, PZ 



(A8) 

(A9) 
-ZP. 



In both cases, T 2 = ±1 defines the ungraded type, i.e. M 
or H. 

To make a correspondence with the AZ scheme, notice 
that T can be labeled as T + or T_, depending on its 
commutation with P. In the even case, this depends on 
a. In the odd case, we have T + = T, while T_ is defined 
as f- Y Z. It follows that f 2 = af 2 . All cases are listed 
in table |H 2 



Appendix B: Semisimple algebras 

For the reader's convenience, we summarize some ba- 
sic facts about semisimple finite-dimensional associative 



2 Note a possibility of confusion: the two methods produce differ- 
ent values of AZ invariants for the same block of a free-fermion 
Hamiltonian because the operator parity is defined by the com- 
mutation with A in one case and H B . P . in the other. This discrep- 
ancy can be fixed by swapping 2+ and T_, or by changing k to 
2 — k. We have chosen not to do that because the variable change 
k i— > 2 — k is convenient in the free-fermion context, see Q|. 
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algebras over R or C. This is quite standard material, see 
e.g. [13 ■ An (associative) algebra A over a field £ is an 
associative ring (i.e., it has addition and multiplication 
with the usual properties) that contains t in its center. 
Here are some examples: 

1. 6(n), the algebra of all n x n matrices whose entries 
are elements of the base field i. It is the same as 
the algebra ?(£) of linear operators acting in the 
space C = t n . 

2. The real numbers (R), complex numbers (C), and 
quaternions (H) regarded as algebras over R. These 
are, in fact, division algebras, i.e. all nonzero ele- 
ments are invertible. According to the Frobenius 
theorem, there are no other real division algebras. 
The only complex division algebra is C. 

3. The group algebra t[G] of a finite group G. It con- 
sists of linear combinations x = ^2 geG x g e g , where 
x g G t, and the basis elements e g satisfy the rela- 
tions e g e h = e gh . 

4. The algebra generated by 1 and e, where e 2 = 0. 
(A general element has the form x + ye, where 
x,y 6 t) 

In this list, items 1 and 2 are simple algebras, 3 is 
semisimple (for 6 = R, C) but generally not simple, and 
4 is not semisimple. 

The simplicity and semisimplicity are defined in terms 
of (finite-dimensional) representations, i.e. vector spaces 
over the base field on which the algebra acts linearly. One 
also needs the notion of intertwinwer: a linear map be- 
tween two representations that commutes with the alge- 
bra action. An isomorphism is an invertible intcrtwiner. 

An algebra A is called semisimple if any subrcpresen- 
tation C C J\f has a complementary representation M., 
i.e. M = C © A4. (This property holds if A is a matrix 
subalgebra that is closed under the Hermitian conjuga- 
tion, one example being the group algebra represented by 
permutation matrices. Indeed, in this case any represen- 
tation Af is unitary, i.e. it has a Hermitian inner product 
such that conjugate elements X, X< G A are represented 
by adjoint operators. Thus, the orthogonal complement 
A4 = C is invariant under the algebra action.) Any rep- 
resentation of a semisimple algebra splits into irreducible 



ones. A semisimple algebra is called simple if it has a 
unique (up to isomorphism) irreducible representation. 

A theorem of Wedderburn states that any semisimple 
algebra is a direct sum of simple algebras, and any simple 
algebra is isomorphic to the algebra of matrices D(n) = 
D®Z(n), where D is a division algebra. D defines a type 
(or equivalence class) of simple algebras. For example, 
the group algebra A = C[G] splits as A = @ a C(£a), 
where C\ arc the irreps of G. In the real case, a similar 
decomposition involves real, quaternionic, and complex 
matrices. 

It is often useful to consider the tensor product of two 
simple algebras over the base field. This is the "multipli- 
cation table" for real division algebras: 

C<g> M C S C©C, (Bl) 
H® M C S C(2), (B2) 
H® R H = R(4). (B3) 

Note that the last two products are simple. That is a 
special case of the following theorem (which also holds 
for base fields other than R and C and can be generalized 
to Z2-graded algebras): If A, B are simple algebras over 
t, and the center of A is equal to t, then A®tB is simple. 

Let A be a simple algebra with center R. While the 
Wedderburn and Frobenius theorems imply that A is iso- 
morphic to a real or quaternionic matrix algebra, there is 
an independent way to characterize A. (This argument 
can actually be used to prove the Frobenius theorem.) 
Consider the complex algebra Ac = A ®r C. It is sim- 
ple, and hence isomorphic to C(n). Thus, the complex 
conjugation o~t ■ Ac — > Ac can be implemented by some 
antilinear operator T , meaning that 

a T (X) = TXf- 1 for all X G .Ac- (B4) 

We may write T = UtK, where Ut is an n x n matrix 
and K is the complex conjugation on C™. It is easy to see 
that T 2 is a nonzero real number. (If the representation 
C is unitary, we may assume that T is anti-unitary so that 
T 2 = ±1.) Now, if f 2 > 0, then A = R(n); otherwise 
A = H(n/2). 
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